Quantum Phase Transitions in Coupled Dimer Compounds 
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We study the critical properties in cubic systems of antiferromagnetically coupled spin dimers 
near magnetic-field induced quantum phase transitions. The quantum critical points in the zero- 
temperature phase diagrams are determined from quantum Monte Carlo simulations. Furthermore, 
scaling properties of the uniform magnetization and the staggered transverse magnetization across 
the quantum phase transition in magnetic fields are calculated. The critical exponents are derived 
from Ginzburg-Landau theory. We find excellent agreement between the quantum Monte Carlo 
simulations and the analytical results. 
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I. INTRODUCTION 



Recent improvements in high magnetic field technol- 
ogy have made detailed investigations of quantum phe- 
nomena in strong magnetic fields possible. In particu- 
lar, they allow for studies of quantum critical proper- 
ties induced by high magnetic fields in weakly coupled 
spin dimer compounds, such as TlCuCl j 1 ! 2 ' 3 , KCuCl^i^, 
BaCuSi 2 Og£, and Sr 2 Cu(B0 3 )2 • • The ground state of 
these materials consists of local spin singlets. If a mag- 
netic field is applied that exceeds their singlet-triplet ex- 
citation gap, they undergo a transition into a magneti- 
cally ordered state. This quantum phase transition into a 
regime with transverse antiferromagnetic (AF) order can 
be described as a Bose-Einstein condensation (BEC) of 
triplet excitations, which behave as bosonic quasiparti- 
cles, called triplons. In the corresponding Bose-Hubbard 
model, this is analogous to the transition from the Mott- 
insulating phase to the superfluid condensate, where the 
magnetic field translates into the chemical potential. At 
very high magnetic fields, there is a saturation thresh- 
old, beyond which the spins fully align along the field 
direction, and the transverse antiferromagnetic order is 
destroyed. A schematic zero-temperature phase diagram 
is shown in Fig^ Because these compounds are three- 
dimensional, field-induced transverse antiferromagnetic 
order persists up to a finite transition temperature, T c Qi), 
between the upper and lower critical field [h c < h < h s ). 
Approaching the critical field from the partially polar- 
ized phase, h — > h c , the critical temperature is expected 
to vanish as T c (h) oc \h — he] 1 /" , with a universal power- 
law exponent that is predicted to be a = 3/2. 7,8 While 
early experimental- and numerical- studies reported non- 
universal, i.e. coupling dependent, values of a > 2, it 
was recently shown in Ref . [l(J that careful fitting of both 
experimental and numerical data indicates an effective 
exponent a{h) > 3/2, which approaches 3/2 as h — > h c . 
Motivated by this observation, Shindo and Tanakaii per- 
formed a fit of specific heat data in TICUCI3 in the close 
vicinity of the lower critical field, yielding an exponent 
of a = 1.67 ± 0.07. This value is closer to the universal 
value of a found in Ref. ^3 than previous experimental 



reports. 

While numerical and experimental results on critical 
properties of the field-induced ordering transition agree 
well with analytical results at finite temperatures&i£, 
there are no comparable scaling predictions for the zero- 
temperature quantum phase transitions. In this pa- 
per, we present a numerical analysis of critical prop- 
erties in cubic systems of coupled dimers. In particu- 
lar, zero-temperature scaling properties of the uniform 
and staggered magnetization in the partially polarized 
phase are studied, based on large-scale quantum Monte 
Carlo (QMC) simulations. These systems, upon increas- 
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FIG. 1: Schematic zero-temperature phase diagram of a 
three-dimensional coupled dimer compound with intra-dimer 
couplings J and inter-dimer couplings J', h denotes the 
strength of the external magnetic field. 

ing the inter-dimer coupling, enter an antiferromagneti- 
cally ordered state at a quantum critical point (J'/J) c . 7 
Magnetic exchange constants are known to depend sen- 
sitively on the distance between the magnetic sites. If 
their relative magnitudes are altered by the application 
of external pressure, such magnetic ordering transitions 
can be induced, as recently observed for TlCuCl3»iii4 In 
this numerical study, the quantum critical point is deter- 
mined for the structures shown in Fig. [21 Furthermore, 
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the quantum criticality induced by an applied magnetic 
field is studied. In particular, the scaling behavior of the 
uniform magnetization and the antiferromagnetic order 
parameter upon entering the partially polarized phase is 
determined. The associated critical exponents can be ac- 
cessed experimentally and thus allow for a quantitative 
comparison. 



H-Hf-H 



fa( fa(' fa( 

M-fai-fat 
fat-H-H 



FIG. 2: Layers of coupled dimers, with staggered (left panel) 
and aligned (right panel) arrangements of dimers. The dimer 
bonds (J) are denoted by solid bars, whereas the inter- 
dimer bonds (J') are denoted by dashed lines. In the three- 
dimensional crystal, these layers are coupled by interlayer 
couplings of the same strength as J'. The arrangement of 
the dimers in neighboring layers is staggered (left panel), or 
aligned (right panel), respectively. 

The paper is organized as follows: in the next section, 
the model used in the numerical simulations is defined 
and details about the quantum Monte Carlo method are 
presented. In Sec. IIIII the zero-temperature phase dia- 
gram of the systems is discussed, and the quantum crit- 
ical points at zero magnetic field are determined. A 
detailed scaling analysis of the magnetic properties at 
the quantum phase transitions in a magnetic field is pre- 
sented in Sec. IIVI Finally, conclusions are given in SeclVl 
The Appendix contains a phenomcnological derivation of 
the mean-field exponents, observed in Sec. IIVI within 
Ginzburg-Landau theory. 



II. MODEL AND METHOD 

We consider the spin-1/2 Heisenberg antiferromagnet 
on the lattice structures shown in Fig. |2 The Hamilto- 
nian is given by 



H 



{id) 



(1) 



where the S, denote localized spin-1/2 moments, and 
Jij indicates the coupling constant between sites i and 
j, which takes values J for the dimer, and J' for the 
inter-dimer couplings. Furthermore, h denotes the ap- 
plied magnetic field. 

Since the lattice structures that are considered here, 
are bipartite (c.f. Fig. [2J, antiferromagnetism in these 
systems is not frustrated, so that their properties can 



be studied using large-scale QMC, without a sign- 
problemi^. Here, we use the stochastic series expan- 
sion QMC metho d 16 i 17 with directed loop update s 18 ! 19 . 
This update scheme results in higher sampling efficiency 
than other methods based on the conventional operator- 
loop update. Close to criticality, autocorrelation times 
are reduced by up to an order of magnitude. This al- 
lows detailed simulations of ground state properties on 
clusters with up to 10 000 sites, even in the presence 
of large magnetic fields. Furthermore, off-diagonal ob- 
servables such as the transverse magnetic structure fac- 
tor can be measured efficiently during the directed loop 
constructionii2i2£ Therefore, the order parameter in the 
partially polarized phase, i.e. the staggered transverse 
magnetization perpendicular to the magnetic field direc- 
tion, can be calculated by using 



(2) 



Here, L denotes the linear system size and the trans- 
verse staggered structure factor, 



(3) 



These observables have been shown to saturate in QMC 
simulations at temperatures below T = J'/2L, which is 
used throughout this study to ensure zero-temperature 
behavior. 



III. ZERO-TEMPERATURE PHASE DIAGRAM 

The zero-temperature phase diagram of coupled dimer 
systems in a magnetic field can be obtained using the 
bond-operator mean-field theory^i^. A schematic phase 
diagram is shown in Fig. 1. It consists of three phases: 
(i) at low fields and small coupling ratios J'/ J, the sys- 
tem is in a magnetically disordered phase, i.e. a dimer 
spin liquid; (ii) at intermediate fields and/or sufficiently 
large values of J'/ J, the ground state is partially spin 
polarized and has an antiferromagnetic long-range order 
transverse to the magnetic field direction; and (iii) at 
large fields, h > h s = J + 5J' all spins are fully po- 
larized. These phase separations occur in both dimer 
arrangements, shown in Fig. |2j 

While bond-operator theory provides a reliable de- 
scription of the phase diagram^, a more precise estimate 
of the critical inter-dimer coupling strength, (J'/J) c , is 
required for the study of critical properties in finite mag- 
netic fields, presented below. In order to determine this 
zero-field quantum critical point, we perform a finite-size 
scaling analysis of the staggered magnetization obtained 
from QMC simulations for various system sizes. Defining 
the dimensionless coupling ratio g = J' /J, the relevant 
finite-size scaling is obtained as follows. The correlation 
length £ diverges near the quantum critical point g c as 
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£ cx \g — 9c\ u ■ The correlation time r c during which 
fluctuations relax and decay (equilibration) , is related to 
the correlation length £ via r c oc £ z cx \g — g c \~" z , with 
the dynamical critical exponent In the vicinity of the 
critical point and at zero-field, the staggered magnetiza- 
tion Trig = m s for g > g c scales as 

m s (x(g-g c ) , (4) 

defining an exponent In general, this implies a finite- 
size scaling relation of the order parameter at the quan- 
tum critical point, g = g c : 

m s cx L-P^. (5) 

Since these systems are explicitly dimerized, the Berry- 
phase contributions to the path-integral cancel, 2 § There- 
fore, this QPT belongs to the universality class of the 
classical 0(3) Heisenberg model in 3+1 dimensions, with 
the dynamical critical exponent z = 1. Since the effec- 
tive classical model is at the upper critical dimension 
(d c = 4), the critical exponents v and r\ take on mean- 
field values v = 1/2, 77 = 0, and the above scaling laws 
hold up to logarithmic corrections^ 

A scaling plot of m^L is shown in Fig. [31 using data 
from system sizes L = 10 to 20. On the scale of the main 
part of Fig. [3] a common intersection point of the finite- 
size data appears to exist. However, on a smaller scale, 
the insets of Fig. [3] exhibit the logarithmic corrections to 
Eq.©. With increasing system size, we observe that the 
intersections for neighboring system sizes move towards 
increasing values of (J'/ J). As seen in inset (b) of Fig. [3 
the systematic increase of the crossing points scales well 
as a function of 1/L 2 , which allows us to extrapolate the 
thermodynamic limit value of (J'/J) c = 0.2492 ± 0.0002 
for the quantum critical point. These logarithmic correc- 
tions are best accessible when system-sizes of different 
orders of magnitude are compared. However, because of 
computational restrictions, we follow the approach de- 
scribed above. 

Performing a similar analysis for the staggered configu- 
ration of dimers in Fig.|^a), we find the quantum critical 
point at (J'/J) c = 0.2170 ± 0.0002. This lower critical 
inter-dimer coupling for the staggered configuration indi- 
cates a reduced tendency towards the formation of dimer 
singlets for this arrangement of dimers, as compared to 
the aligned dimer configuration. 

Beyond the critical inter-dimer coupling, the staggered 
magnetization increases with a mean-field exponent (3 = 
1/2, as in Eq. gj). The QMC data of m s , shown in Fig.H 
are consistent with this scaling behavior. However, due 
to restricted system sizes, we are not able to extract the 
logarithmic corrections to the scaling law. 

IV. FIELD-INDUCED QUANTUM PHASE 
TRANSITION 

Having determined the quantum critical points of 
coupled-dimer arrays, we proceed to study the effects of a 
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FIG. 3: Scaling plot of the zero-temperature staggered mag- 
netization in the aligned arrangement of Fig. 2 (b) obtained 
from quantum Monte Carlo simulations using systems of lin- 
ear sizes L = 10 to 20. J denotes the intra-dimer coupling. At 
the critical inter-dimer coupling (J'/J) c , the different curves 
intersect each other (main panel). Inset (a) shows a mag- 
nification of the intersection region. Inset (b) exhibits the 
corrections to the mean-field scaling behavior, arising from 
logarithmic corrections in the upper critical dimension. The 
statistical error bars fall within the symbol size. 
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FIG. 4: Staggered magnetization m s of coupled dimer arrays 
as a function of the inter-dimer coupling strength J', for the 
two different dimer arrangements of Fig. 2. Dashed lines 
show the mean-field scaling behavior with exponent f) = | 
close to the quantum critical point. 



magnetic field in the different regimes. In particular, we 
determine the scaling behavior of the uniform magnetiza- 
tion and the order parameter upon entering the partially 
polarized region. The corresponding scaling exponents 
are accessible experimentally by direct measurement of 
the magnetization, and by neutron scattering. Thus, the 
numerical results can be compared with analytical pre- 
dictions and with measurements on the materials men- 
tioned above. Here, simulations for the aligned configu- 
ration of dimers, shown in Fig.^b), are presented. For 
the staggered configuration of dimers, the same scaling 
exponents are obtained. This is due to the underlying 
universality. 
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A. Scaling of the Uniform Magnetization 

First, we discuss the behavior of the uniform magneti- 
zation m u (h) as a function of the applied magnetic field. 
In the dimer spin liquid phase, i.e. for coupling ratios 
smaller than ( J'/ J) c , a finite magnetic excitation gap A 
separates the ground state singlet and the lowest triplet 
state. Thus, a finite magnetic field h c — A is required to 
close this gap, and to induce a finite uniform magnetiza- 
tion. The Ginzburg-Landau approach presented in Ap- 
pendix A predicts m u to increase linearly, m u oc (h — h c ) 
for h > h c = A and for J' < J' c . When the inter-dimer 
coupling is increased, the gap A becomes smaller, until 
it vanishes at the critical inter-dimer coupling J' c . At the 
quantum critical point, J' = J' c , the Ginzburg-Landau 
approach predicts the uniform magnetization to scale as 
m u oc h 3 . For larger values of J' > J' cl the excitation 
gap remains zero, and the finite uniform susceptibility 
Xu results in a linear response m u = Xuh of the uniform 
magnetization in the Neel-ordered regime beyond J' c . 

Let us now compare these predictions with the QMC 
results obtained from simulations of systems with linear 
sizes up to L = 30. For such large system sizes, we did 
not detect finite-size effects in the uniform magnetiza- 
tion. In Fig. [31 results for various values of the inter- 
dimer coupling strength are presented. The main part 
of Fig. UJ shows the uniform magnetization over the full 
range of magnetic field strengths, up to the saturation 
field h s = J + 5 J'. Inset (a) of Fig. [3] focuses on the 
region close to the critical field h c for various inter-dimer 
couplings J'. The linear scaling in h, which is predicted 
on both sides of the quantum critical point, is clearly 
observed in this small-m^ region. In contrast, at the 
quantum critical point, J' = J' c , the uniform magnetiza- 
tion increases non-linearly with h, as shown in Fig.[5J(a). 
Indeed, we observe a scaling m u oc h 3 , presented in 
Fig. UJ (b) . This is expected from the Ginzburg-Landau 
theory (c.f. Appendix A) and the bond-operator mean- 
field theory— Using QMC, these scaling exponents of 
the uniform magnetization are verified. 

Finally, we note that the field dependence of m u for 
the entire region between h c and h s is linear for weakly 
coupled dimers (see e.g. Fig.[S]for J' = 0.07 J), consistent 
with bond-operator mean-field theory^. For larger val- 
ues of J', m u shows deviations from this linear behavior 
in high magnetic fields. This non-linear behavior can be 
accounted for within bond-operator theory by including 
the contributions of higher-energy triplet modes to the 
ground state—. 

B. Scaling of the Order Parameter 

Next, we discuss the scaling properties of the order 
parameter in the partially polarized phase, i.e. the 
staggered transverse magnetization perpendicular to the 
magnetic field direction, mjr(h), as a function of the ap- 
plied magnetic field h. Before presenting numerical data, 
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FIG. 5: Zero-temperature uniform magnetization of aligned 
dimer arrays for different inter-dimer couplings as a function 
of the magnetic field h. Inset (a) focuses on the low-field 
region, h/J < 1. The scaling m u oc h 3 at J' = J' c is demon- 
strated in inset (b). 



the expectations from the Ginzburg-Landau theory from 
Appendix A can be summarized as following. For J' < J' c 
and magnetic fields h > h c = A, mj" oc \J (h — h c ) is 
expected, consistent with bond-operator mean-field the- 
ory.— At the critical inter-dimer coupling, J' = J' c , a 
linear relation mjr oc h is expected. Within the antiferro- 
magnetically ordered regime, J' > J' c , a finite staggered 
magnetization mjr = m s exists. For small fields, the or- 
der parameter scales with mjr — m s oc h 2 , as presented 
in Ref . |H 

QMC calculations of mjr for different values of J'/ J in 
a system with L — 16 are shown in Fig. This Figure 
demonstrates that for J' < J' c , the staggered magnetiza- 
tion is largest for h = (h c + h s )/2 and decreases upon 
approaching h c and h s . In the thermodynamic limit, mjr 
vanishes in the fully polarized regime, h > h s and below 
h c for weakly coupled dimers. In contrast to the uni- 
form magnetization, the field dependence of shows a 
strong finite-size dependence. Finite-size effects are most 
pronounced for magnetic field regimes outside the par- 
tially polarized phase, i.e. for h < h c and h > h s , where 
mjr = in the thermodynamic limit. On the finite sys- 
tems accessible to QMC, the values of nig , defined in 
Eq. (0) do not vanish. However, they scale to zero upon 
increasing the system size. In order to extract the scal- 
ing behavior of , we thus need to perform a careful 
finite-size scaling analysis of the numerical data. 

Such an analysis can be directly applied in the weak 
coupling regime, for J' < J' c . Here, the field- induced or- 
dering transition is known to constitute a Bosc-Einstein 
condensation of triplonsji with a dynamical critical expo- 
nent z = 2, reflecting the quadratic dispersion relation of 
these bosonic excitations. Therefore, the quantum phase 
transition at the critical field h c is in the universality 
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FIG. 6: Zero-temperature staggered transverse magnetiza- 
tion mj~ in coupled aligned dimer arrays as a function of an 
applied magnetic field h. Quantum Monte Carlo data of a 
16 x 16 x 16 system are shown for different values of the inter- 
dimer coupling J'. The intra-dimer coupling is denoted by J. 
The off-sets of raj: outside the range h c < h < h s are due to 
finite-size effects. 




FIG. 7: Scaling plot (a) and data collapse (b) for the 
zero-temperature staggered transverse magnetization mj: of 
weakly coupled aligned dimers as a function of the applied 
magnetic field h for a coupling-ratio of J' /J = 0.2. Results 
from quantum Monte Carlo simulations are shown for systems 
sizes of L = 10 to 20. 



class of the classical 5-dimensional 0(2) model, and a 
finite-size scaling analysis similar to the one performed 
for the zero-field quantum critical point in Sec. IIV Al 
can be applied. Since the classical theory is now above 
d c = 4, mean-field scaling without logarithmic correc- 
tions applies. In particular, at the critical field /i c , one 
finds 



L cx F(L/0 



(6) 



where F is a scaling function that depends on the ratio of 
the correlation length £ and the system size L. In order 
to determine h c for a given value of J' < J c , we can thus 
construct a scaling plot of m^L as a function of h/J. 
Fig. illustrates an example of such a plot for an inter- 
dimer coupling J' — 0.2J. In Fig. [7{a), data for rajrL 
are shown for system sizes L = 10, 14, 18, and 20. The 
intersection point of the different finite-size data allows 
the extraction of the critical field h c = 0.383 ± 0.001 for 
this value of J' = 0.2 J. 

Once the critical field h c is determined, data collapse 
of the finite-size data is verified in the vicinity of h c . 
Namely, from the scaling of the correlation length close 
to the critical field, £ cx \h — h c \~", one obtains 



L cx F(L 1/lJ \h-h c \ ) 



(7) 



polynomials. Using these extrapolated values of mjr(h), 
representing the thermodynamic limit, we then obtained 
the scaling behavior of m^{h) shown for different values 
of J'/ J in Fig. ED 
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with a new scaling function F and v = 1/2. An example 
of the data collapse for J' = 0.2 J is shown in Fig. 0(b). 

In order to extract the behavior of mjr(h) for h > h c 
in the thermodynamic limit, we perform an extrapolation 

of the finite-size data mjr{L) as a function of 1/L. An nent | is also found for stronger inter-dimer couplings, 



FIG. 8: Extrapolation of finite-size data of staggered mag- 
netization for different magnetic fields close to the transi- 
tion. At h c , thermodynamic limit is obtained from a 1/L- 
extrapolation. Higher order polynomials are required away 
from criticality. For all data presented here, up to the third 
order is used. Dashed lines are for h < h c . 

For weakly coupled dimers with J' = 0.07 J Fig. EJa) 
exhibits a scaling of cx \/h — h c . The scaling expo- 



example of such an extrapolation, again for J' / J = 0.2 
is shown in Fig. [S] for various values of h. While for h 
close to h c , a linear scaling in 1/L is obtained, extrap- 
olation away from h c requires the use of higher order 



as shown in Fig.|5fb) for J' = 0.2 J. This square root de- 
pendence cx \Jh— h c is consistent with the analytical 
results and expected from the mean- field value f3 = 1/2. 
Next, let us consider the scaling of m^{h) for a criti- 
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cal inter-dimer coupling J' = J' c . In Fig, [ffic) . we show 
results from simulations at the critical point J' = 0.249. 
Here, a linear scaling oc h is observed, consistent with 
the Ginzburg-Landau calculations from Appendix A. Fi- 
nally, we consider the Neel-ordered phase (J' > J^). The 
results from the extrapolated QMC data is shown for 
J' = 0.3 J in Fig. EJd). Increasing the magnetic field, 
a non-linear increase of mjr(h) from its zero-field value 
m s is found, which can be fitted well to the analytical 
prediction, mjr(h) — m s oc h 2 . We thus obtain agreement 
of the QMC data with the analytical results based on 
Ginzburg-Landau and bond-operator mean-field theory 
also for the scaling of the order parameter. 
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FIG. 9: Scaling behavior of the zero-temperature staggered 
transverse magnetization as a function of the applied mag- 
netic field h. Quantum Monte Carlo data are shown after 
extrapolation of finite-size data (L = 10 — 20) to the thermo- 
dynamic limit. (a)(b) in the dimer spin liquid phase (J' < J' c ), 
(m^) 2 scales linear with h — h c . (c) near the critical coupling 
(Jc), m s x scales linear with h — h c . (d) for J' > J' D , — m a 
increases quadratically with h. 



V. SUMMARY 

We examined quantum phase transitions in three- 
dimensional coupled dimer arrays. The zero-temperature 
phase diagrams of these systems feature a low-field dimer 
spin liquid phase at weak inter-dimer couplings, a par- 
tially polarized regime with long-range transverse mag- 
netic order for intermediate magnetic fields h c < h < h s , 
and a fully polarized phase at high magnetic fields. 

The critical exponents associated with quantum phase 
transitions between these regimes were extracted us- 
ing finite-size scaling analysis of quantum Monte Carlo 
data. The numerical values of these exponents com- 
pare well with Ginzburg-Landau calculations and bond- 
operator mean-field theory. In particular, for small inter- 
dimer coupling, the order parameter is found to scale as 
{h — hc) 1 / 2 , whereas the uniform magnetization scales lin- 
ear in h — h c . Moreover, at the quantum critical point 
(J'/J)c w e demonstrated linear scaling for mjr oc h. In 
a recent magnetization study on the TICUCI3 system un- 
der hydrostatic pressure^, a cubic scaling of the uniform 
magnetization m u oc h 3 at the critical value of the ap- 
plied pressure was observed. This is in perfect agreement 
with numerical results and analytical predictions based 
on Ginzburg-Landau calculations and bond-operator the- 
ory. 

We found that corrections to mean-field scaling emerge 
at the zero-field pressure- induce d quantum phase transi- 
tion, such as those in Ref. 113114 For example, in Ref. IhH 
the pressure-dependence of the gap was fitted to a power 
law with an exponent 0.445, while the mean-field theory 
predicts an exponent 0.5, without taking logarithmic cor- 
rections into account^ Assuming a linear scaling of the 
pressure P with the inter-dimer interactions^ close to 
the critical point P c , we find the scaling of the gap A to 
be 



Acx 1/f oc {P - P c 



a/2 



In 



P-Pr 



P, 



1/6 



(8) 



from a renormalization group analysis of the classical <fi A - 
theory at d c = 4 26 . Therefore, careful fitting of experi- 
mental data can exhibit the logarithmic corrections near 
quantum critical points in three-dimensional compounds. 
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APPENDIX A: GINZBURG-LANDAU THEORY 

In this Appendix, we present a phenomenological 
Ginzburg-Landau description of weakly coupled dimers 
in a magnetic field 28 , and derive the scaling exponents of 
the uniform magnetization m u and the order parameter 
mjr in the partially polarized phase. 

Due to rotational symmetry within the plane perpen- 
dicular to the external field, the free energy of the system 
can be expanded as 



F = a(h,J') (m^) 2 + b (m^) 4 , 



(Al) 



neglecting higher orders of mj" in the vicinity of the 
critical fields Here, a(h, J') vanishes for h equal to 
the critical magnetic field h = h c (J'). Taking into ac- 
count time reversal symmetry, a can be expressed as 
a = a ■ (h 2 — /i 2 (J')). The staggered magnetization 
and the uniform magnetization m u is obtained from 
Eq. Ej]as 



OF 



dm 4 

h,J> 

From this, one finds 



and m u 



dF 
~dh 



(A2) 



To proceed further, we distinguish two regions, namely 
weak inter-dimer couplings (J' < J' c ) and the quantum 
critical point (J' = J' c ). For the first region, (|A3|) can be 
expanded in terms of the distance to the lower critical 
field, i.e. h — h c . One obtains 



for J' < J' 



a(h + h c )(h — h c ) 



2b 

a/2 



oc (h — h c 
■ m u tx 2ah(h — he) 1 



(A5) 
(A6) 



At the critical inter-dimer coupling the critical field is 
zero. Therefore, 



for J' = Jt 



' ' 2b 
m„ oc 2ah 3 . 



(A7) 
(A8) 



m a = 



~a(h? - hi) 



2b 



2ah(rrig 



(A3) These exponents agree with the numerical data pre- 

sented in the preceding sections, as well as previous bond- 
(A4) operator mean-field theory results 2 ^. 
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